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1  EXECUTIVE  SUMMARY 

The  ultimate  goal  of  proposed  research  is  to  develop  a  new  advanced  simulation  and  designing 
tool  of  textiles  and  textile-based  composites  for  developing  lighter,  more  durable  and  affordable 
composite  structures  in  military  and  commercial  applications.  The  mechanism-based  damage 
model  developed  in  this  research  enables  predicting  validated  physics  of  failure  characterization 
of  textile-based  composites  so  that  eliminates  conservative  design  penalties.  Different  from  the 
existent  composite  models,  the  new  model  addresses  issues  related  to  different  local  damage 
mechanisms  based  on  “meso-structural”  behaviors  of  fibers  with  the  use  of  a  truly  nonlinear  fiber 
model,  high-fidelity  fiber-to-fiber,  fiber-to-other  material  interaction  models,  and  an  enhanced 
homogenization  method.  The  developed  model  can  greatly  enhance  accuracy  and  efficiency  of 
simulation  and  design  for  textile-based  composite  structures,  such  as  for  advanced  composite 
armors,  under  various  static  and  dynamic  loading  conditions,  while  further  extension  can  also 
include  ballistic  and  blast  loading  conditions.  The  developed  innovative  textile  mechanics  model 
and  simulation  and  designing  tool  will  also  find  important  applications  in  protective  clothing  and 
parachute  applications  for  protecting  individual  soldiers. 

In  the  current  program  of  the  proposed  research,  we  have  developed  and  integrated  the 
following  two  new  capabilities  to  our  simulation  and  design  system: 

1.  Extended  nonlinear  fiber  model  for  simulating  general  textile  fabrics  and  textile 
composite  materials. 

2.  A  high  fidelity  meso-structural  model  for  textile  composites. 

We  have  extended  our  nonlinear  fiber  model  to  general  three-dimensional  fabric  system  with 
arbitrary  configuration  and  with  interactions  with  other  fibers  and  materials  in  the  composite  for 
simulating  a  practical  textile  composite.  The  new  model  accounts  for  geometry  nonlinearity  and 
arbitrary  large  deformation  of  fibers,  fiber-to-fiber  contacts,  and  fiber-to-other  material  contacts 
in  a  textile  composite.  This  results  in  a  unique  capability  for  simulating  textile  composites, 
especially  for  more  accurate  damage  assessment.  The  new  capability  will  be  further  extended  and 
integrated  with  the  other  advanced  capabilities  being  developed  for  designing  truly  optimum 
textile  composites  for  military  and  commercial  applications. 

Keywords: 

Textile  Composites,  Lightweight  Structure  Technologies,  Future  Combat  Systems,  Tactical  and 
Support  Vehicle  Fleets,  Mechanism-based  Damage  Model,  Textile  Fabrics. 
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2  BACKGROUND 

One  of  Army  research  goals  within  the  next  ten  years  is  to  develop  design  and  analytical  tools 
that  can  predict  the  impact  damage  and  the  corresponding  post-impact  performance  and  strength 
of  the  composite  structures.  Another  special  interest  the  Army  has  is  to  develop  new  experimental 
and  analytical  methodologies  for  textile  materials  that  are  used  for  protective  clothing  and 
parachutes  for  the  soldiers.  The  objective  of  proposed  research  is  to  develop  advanced  simulation 
and  designing  tools  for  textiles  and  textile-based  composites  for  the  need  to  design  lighter,  more 
durable  and  affordable  composite  structures  for  the  Future  Combat  Systems  (FCS)  and  other 
military  applications.  The  mechanism-based  damage  model  proposed  in  this  research  will  enable 
validated  physics  of  failure  characterization  to  eliminate  conservative  design  penalties.  Different 
from  existing  models,  the  new  model  developed  in  this  program  will  address  issues  of  different 
local  damage  mechanisms  based  on  “meso-structural”  behaviors  of  fibers,  truly  nonlinear  fiber 
model,  high-fidelity  fiber-matrix  interaction  model,  and  an  enhanced  homogenization  method. 
The  developed  model  can  greatly  enhance  simulation  and  design  of  the  textile  composite 
structures  in  ground  vehicles  under  various  loading  conditions,  including  ballistic  and  blast 
impacts.  The  innovative  textile  mechanics  model  and  simulation  and  designing  tools  will  also 
find  important  applications  in  protective  clothing  and  parachutes  for  the  soldiers. 

We  are  accumulating  advanced  and  unique  analytical  and  design  capabilities  for  developing 
innovative  composite  structures  for  Army  future  combat  and  tactical  vehicles.  The  new 
capabilities,  developed  and  being  developed,  include: 

1.  Function-oriented  material  design  (FOMD)  with  a  multi-domain,  multi-step,  multi¬ 
material  topology  optimization  process. 

2.  Multi-scale  modeling  and  homogenization  method. 

3.  Fabrics  mechanics-based  damage  prediction. 

4.  Truly  nonlinear  fiber  model,  which  accounts  for  arbitrary  large-deformation  of  elastic 
fiber,  uneven  tension  distribution  along  the  fiber,  fiber-fiber  or  fiber-matrix  interaction 
for  high-fidelity  prediction  of  fiber  behavior  in  a  composite  material. 

5.  Fiber  arrangement  optimization  in  composite  materials. 

6.  Advanced  structural  analysis  capabilities  for  durability,  NVH,  crashworthiness,  ballistic 
and  blast  simulations. 
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These  advanced  capabilities  are  being  integrated  into  a  unique  simulation  and  design  platform, 
which  can  be  used  for  designing  and  evaluating  innovative  composite  materials  and  structures  in 
future  vehicle  systems.  Following  are  brief  summaries  of  the  methodologies  being  developed. 

2.1  Advanced  Textile  Composite  Damage  Modeling 

Textile-based  composites  are  outstanding  candidates  for  lightweight  vehicle  structural  materials. 
Textile  fabrics  are  also  widely  used  in  soldier  systems,  such  as  protective  clothing,  parachutes, 
and  shelters.  Predictions  of  elastic  properties  of  textile  composites  can  now  be  found  in  many 
literatures.  However,  the  non-linear  behaviors  and  the  damage  mechanism  of  textiles  and  textile- 
based  composites  have  not  been  fully  understood  due  to  the  complex  geometry  effects  and 
anisotropic  non-linear  properties  of  the  associated  meso-structures.  As  a  result,  damage  modeling 
of  textile  composites  is  still  an  active  research  area.  There  are  various  models  for  textile 
composite  damage  simulation,  which  either  greatly  simplify  the  material  behaviors,  or  are  too 
complicated  and  inefficient  in  aspect  of  practical  usage.  In  the  previous  models,  the  textile 
composites  are  modeled  as  materials  similar  to  homogeneous  materials.  Especially  in  damage 
modeling,  the  detailed  damage  sequences  and  mechanisms  are  not  actually  accounted  for.  In  our 
proposed  model,  textile  composites  will  have  detailed  texture  and  “meso-structures”  with  the 
consideration  of  essential  fiber  nonlinear  characteristics,  which  can  be  used  to  predict  the  meso- 
structural  behaviors  of  the  textile  and  textile  composites,  including  delaminations  and  local 
damages.  An  enhanced  homogenization  method  will  be  then  developed  to  include  this  textile 
composite  model  for  the  macro-structural  analysis  with  the  needed  computational  efficiency.  This 
results  in  a  multi-scale  modeling  process,  which  can  be  used  to  predict  both  “global”  and  “local” 
behaviors  of  the  composite.  The  major  focus  is  to  develop  an  advanced  fiber-matrix  model  for  the 
textile  composite.  This  model  will  base  on  the  local  mechanical  behavior  of  the  material  and  local 
damage  mechanism,  so  that  high  fidelity  damage  prediction  can  be  obtained.  The  proposed  multi¬ 
scale  model  will  also  be  highly  efficient  for  simulation  and  design  of  composite  due  to  the 
enhanced  homogenization  method  and  advanced  computational  algorithms  to  be  developed. 
Figure  1  summarizes  the  scope  of  work  for  the  research  at  The  University  of  Michigan  related  to 
textile  composites. 
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Problems  Involved  and  Modeling  Work: 


Figure  1.  Scope  of  work  of  the  research  at  U  of  M  related  to  textile  composites 


2.2  Simulation  and  Design  Capabilities  to  be  Integrated 
2.2.1  Function-oriented  material  design  (FOMD) 

We  have  developed  an  award-winning  advanced  CAE  tool,  called  FOMD  (Function-Oriented 
Material  Design)  [1-14],  which  employs  a  multi-domain,  multi-step  and  multi-material  topology 
optimization  approach.  Our  FOMD  technology  has  received  a  best  paper  award  in  2003  from  the 
ASME  Design  Engineering  Division  [1].  FOMD  has  been  developed  as  a  commercial-grade 
CAE  tool  for  designing  advanced  and  high-performance  structures  and  materials.  Figure  2 
illustrates  example  graphic  user  interfaces  of  a  FOMD-Blast  simulation  and  design  system,  which 
is  being  developed  for  designing  advanced  blast-protective  composites  [14].  As  shown  in  Fig.  2, 
this  software  system  includes:  1)  explosive  modeling,  2)  soil  modeling  for  landmine-soil 
interaction,  3)  geometry  modeling  of  vehicle  structure  and  occupants,  4)  blast  prediction,  which 
will  predict  both  blast  pressure  and  fragment  penetrations  on  the  composite,  5)  composite  layout, 
which  includes  raw  material  selection  and  basic  material  configuration,  6)  optimal  design  based 
on  the  FOMD  and  topology  optimization.  This  software  system  will  be  continuously  extended  to 
include  more  advanced  modeling  and  simulation  capabilities  and  post-processes  that  will  finalize 
the  designs  for  prototyping  and  mass  production. 


5 


Explosive  modeling 


Optimal  design 


Soii  modeling 


Geometry  modeling 


Composite  layout 


Blast  prediction 


Figure  2.  Proposed  simulation  and  design  process 


2.2.2  Multi-scale  modeling  and  homogenization  method 

Figure  3  illustrates  fundamental  idea  of  the  multi-scale  modeling.  As  shown  in  Fig.  3,  a 
completed  modeling  process  should  starts  from  the  micro-level  with  the  consideration  of  micro- 
structural  behavior  of  the  composite  material.  For  example,  detailed  micro-structure  of  a  fiber  or 
thread  can  be  considered  at  this  stage.  Then  effective  material  properties  of  the  fiber  or  thread  at 
the  meso-level  can  be  estimated  through  a  homogenization  process  from  the  micro-level  to  the 
meso-level.  Next,  a  higher  level  homogenization  process  can  be  employed,  which  estimates 
effective  material  properties  of,  for  example,  a  cut-off  piece  of  the  composite,  from  the  meso- 
level  to  the  macro-level.  Finally,  a  structural  analysis  can  be  conducted  at  the  macro-level 
without  considering  composite  details  at  the  micro-  and  meso-  levels.  This  micro-meso-macro 
modeling  process  makes  the  composite  analysis  highly  efficient  and  effective  without  loosing  any 
detail  of  the  composite  at  the  meso-  and  micro  levels.  In  next  step,  we  will  develop  and 
implement  a  new  generated  homogenization  method,  which  will  suit  the  needs  described  here. 
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Non-Homogenized  Material  Homogenized  Orthotropic  Material  Fiber  Bundle 


Fiber  Bundle  Braided  Structure  Homogenized  Material 


Figure  3.  Hierarchical  modeling  at  micro-,  meso-,  and  macro-levels 
2.2.3  Fabrics  mechanics-based  damage  prediction 

An  advanced  model  for  fabrics  has  been  developed  at  the  University  of  Michigan  with  the 
assumption  of  small  fiber  deformation  and  linear  geometry  [15].  This  model  is  unique  compared 
to  the  state-of-art  models  in  several  aspects:  1)  it  is  based  on  physical  behaviors  of  fiber  yams 
instead  of  assumptions;  2)  it  is  numerically  efficient;  3)  it  enables  automatic  modeling  for  textile 
composites  with  a  few  design  parameters;  4)  it  can  be  modified  to  use  in  textile  modeling.  Figure 
4  illustrates  some  example  geometrical  modeling  works  conducted  at  the  U  of  M  with  a 
traditional  finite  element  approach  [16]. 
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Figure  4.  Fabrics  mechanics  and  geometrical  modeling  [2] 

There  are  different  mechanisms  in  textile  composites  damage  process.  For  example,  the 
damage  of  woven  composites  will  include  the  delamination  of  fiber  tows,  and  the  sliding  between 
the  separated  tows  after  the  initial  damage.  Most  current  models  cannot  address  these  damage 
mechanisms.  In  the  framework  of  our  general  nonlinear  multi-scale  model,  this  kind  of 
mechanism  can  be  easily  considered,  and  the  influences  of  different  damage  mechanism  on 
overall  material  performances  can  be  investigated  in  details.  The  freedom  provided  will  also 
include  the  possibility  to  consider  additional  new  damage  mechanisms  observed  in  testing 
processes. 

3  TECHNICAL  OBJECTIVES  AND  PROPOSED  TASKS 
3. 1  Technical  Objectives  of  the  Overall  Research 

The  ultimate  goal  of  the  proposed  research  is  to  develop  a  comprehensive  and  advanced 
simulation  and  designing  tool  for  textiles  and  textile-based  composites  for  the  need  of  designing 
lighter,  more  durable  and  affordable  composite  structures  for  future  military  applications.  The 
mechanism-based  damage  model  proposed  in  this  research  will  enable  validated  physics  of  failure 
characterization  to  eliminate  conservative  design  penalties.  The  new  model  developed  in  this 
program  will  address  issues  of  different  local  damage  mechanisms  based  on  “meso-structural” 
behaviors  of  the  fibers.  We  are  proposing  to  develop  and  integrate  the  following  new  capabilities 
to  our  simulation  and  design  system: 
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1 .  Innovative  mechanical  model  for  textile  fabrics. 

2.  Nonlinear  multi-scale  modeling  for  textile  composites. 

3.  Mechanism-based  damage  modeling  for  high-fidelity  prediction. 

Detailed  technical  objectives  of  the  overall  research  are  enumerated  as  follows: 

1 .  Extend  our  nonlinear  fiber  model  for  simulating  general  textile  fabrics  and  composite 
materials 

2.  Develop  high  fidelity  meso-structural  model  of  textile  composites 

3.  Develop  advanced  homogenization  theory  for  micro-meso-macro  modeling  of  textile 
composites 

4.  Develop  advanced  mechanism-based  damage  theory  for  textile  composites 

5.  Develop  effective  and  efficient  numerical  tools  for  analysis  and  design  of  textile 
composites 

6.  Conduct  proof-of-concept  through  prototyping  and  example  applications 

In  the  current  program,  we  proposed  to  focus  on  the  first  two  tasks  as  a  preliminary  and 
feasibility  study  for  the  further  developments  in  this  research. 

3.2  Proposed  Tasks  in  the  Current  Program 

The  following  tasks  were  proposed  to  be  carried  out  in  the  current  phase  of  the  research  program: 

3.2.1  Task  I  -  Extend  our  nonlinear  fiber  model  for  simulating  general  textile 
fabrics  and  textile  composite  materials 

In  this  task,  we  proposed  to  extend  our  nonlinear  fiber  model  to  general  three-dimensional  fabric 
system  with  various  configurations  for  simulating  practical  textile  fabrics  and  composite 
materials.  The  formulation  will  account  for  geometry  nonlinearity  and  arbitrary  large  deformation 
of  the  textile  fabrics.  This  will  result  in  a  unique  capability  for  simulating  textile  composite 
materials,  especially  for  the  damage  assessment. 

3.2.2  Task  II  -  Develop  a  high  fidelity  meso-structural  model  for  textile 
composites 

In  this  task,  we  proposed  to  develop  a  high-fidelity  meso-structural  model  for  general  textile 
composite,  which  will  include  fiber/fabrics  and  filling/matrix  material.  This  meso-structural 
model  would  allow  us  to  conduct  a  detailed  localized  analysis  for  the  meso-behaviors  of  the 
composite,  which  include  nonlinear  deformation  and  strength  of  individual  fibers  in  the 
composite. 
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4  MAJOR  RESULTS 

4. 1  Development  of  the  New  Nonlinear  Fiber  Model 


We  have  developed  a  new  capability  to  simulate  a  truly  nonlinear  and  large  deformation  fiber 
system  based  on  an  efficient  and  advanced  formulation  of  three-dimensional  fiber  systems.  [15] 
The  new  fiber  model  is  in  a  very  compact  form,  resulted  from  employing  proper  coordinate 
systems  and  state  variables.  It  accounts  for  geometry  nonlinearity,  extensibility  of  the  fiber, 
tension  variation  along  the  fiber,  tension-dependent  mechanical  properties  of  the  fiber,  and  other 
nonlinear  effects  in  such  fiber  system  undergoing  arbitrary  large  deformation. 

4.1.1  Basic  equations 

4. 1.1.1  State  equations 

Figure  5 -a  illustrates  a  general  thread  member  taken  out  from  a  fabric  system  with  the 
adoption  of  a  curvilinear  coordinate  system.  Figure  5-b  illustrates  the  equilibrium  of  an 
infinitesimal  element  cut  out  from  the  thread  at  point  P.  As  shown  in  Fig.  5-a,  for  an  arbitrary 

material  point  P  on  the  thread,  vector  r  =  |x,  y,  z}^  represents  the  position  of  P  in  the  global 

coordinate  system,  O-XYZ.  In  general,  r  is  a  function  of  space  variable  and  time  variable.  For 
deriving  the  thread  model  with  an  arbitrary  initial  shape  and  potential  large  deformation,  we 

adopt  a  local  coordinate  s,  which  represents  undeformed  arc-length  oP  measured  from  the  left 
end  of  the  thread  (namely  point  o  in  Fig.  5)  to  the  material  point  P  along  the  axial  curve  of  the 


thread. 


Y 


Y 


z 


a)  A  thread  member 


b)  Equilibrium  of  an  infinitesimal  element 


Figure  5.  Modeling  of  a  nonlinear  thread  member 
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Figure  5-b  illustrates  the  equilibrium  of  the  infinitesimal  element  cut  out  from  the  thread  at 
point  P.  Here,  ds  denotes  the  undeformed  arc  length  of  the  infinitesimal  element,  d p  denotes 
the  elongated  arc  length  of  the  same  element  after  fiber  deformation  due  to  the  tension  and 

external  force  applied  on  the  element,  T  =  T  |  denotes  the  tension  vector  acting  on  the 

point  P,  f  =  {/x?  Z;;? /z}  denotes  the  external  force/per  length  acting  on  the  element,  and 


m — ^  is  the  inertial  force/  per  length,  where  m  denotes  density  of  the  thread,  and  a  =  — ^  is 


dt 


dt 


the  acceleration  vector  of  the  thread  element.  It  is  assumed  that  the  external  force  is  non¬ 
following  with  a  fixed  coordinate  5  at  the  arc-length  coordinate  system,  and  its  amplitude  and 
direction  are  independent  of  the  elongation  of  the  thread,  namely  we  have  f  =  .  Note  that  a 


following-type  force  can  also  be  considered  and  corresponding  formulation  can  be  obtained  using 
a  similar  approach  described  in  this  report. 


From  Fig.  5-b,  the  equilibrium  of  the  forces  acting  on  the  thread  element  can  be  obtained  as 


5T  aV  „  . 

- m — ^  +  f  =  0  . 

ds  dt^ 


(1) 


Assuming  that  T  is  the  magnitude  of  the  tension  vector  T,  and  n  is  the  unit  direction  vector 
of  T,  we  have, 


T  =  7h  (2) 

where  T  can  be  calculated  using 

r=A^.  ,3, 

From  Fig.  5-b,  the  direction  vector,  n,  can  be  calculated  as 


(4) 


here,  elongation  of  thread  element  has  been  considered  in  Eq.  (4),  this  results  in  a  formulation 
that  can  account  for  extensibility  of  the  fiber  or  thread. 


To  simplify  the  problem,  let  us  consider  that  the  thread  member  is  made  of  a  linear  elastic 
material,  which  obeys  a  constitution  law: 
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T  =  EA 


(5) 


yds  y 

which  implies  that  the  tension  force  T  is  determined  by  tensile  modulus,  E,  thread  cross-section 
area,  A,  and  strain  (elongation/length),  namely  {dp- ds^jds .  From  Eq.  (5),  we  can  obtain 


ds  _  EA 
d  p  EA  -l-  T 

Substituting  Eq.  (4)  into  Eq.  (2)  and  using  derivative  chain  rule  obtains 

dp  ds  dp 

Substituting  Eq.  (6)  into  Eq.  (7)  obtains 

T=f  — 
ds  ’ 

where 


(6) 


(7) 


(8) 


^  EA-T 

T  = - 

EA  +  T 


(9) 


is  called  generalized  tension  in  this  report,  which  has  the  same  dimension  as  the  tension.  Note  that 
for  an  inextensible  thread,  we  have  EA  oo  and  therefore  T  ^  T  . 


Substituting  Eq.  (8)  into  Eq.  (1),  the  governing  equation  for  the  thread  member  can  be 
obtained  as 


ds 


y  ds 


d\  „ 
+  m — ^  =  t  . 

de 


(10) 


It  is  seen  that  the  equation  of  motion  becomes  very  simple,  which  is  resulted  from  the 
introduction  of  a  new  generalized  tension  variable,  f ,  and  the  local  arc-length  coordinate. 


For  a  static  equilibrium  problem,  assuming  that  f  =  is  static  force  acting  on  the  thread, 
Eq.  (10)  can  be  reduced  as 
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(11) 


ds 


,dY 

ds 


=  ^s, 


where  T  =  r(r)  is  a  function  of  the  thread  configuration  through  the  relationship  defined  by  Eqs. 


(9),  (5),  and  ^=1^ 

CIS  V  V 


ds 


Vibration  problem  of  the  thread  is  usually  considered  about  the  static  equilibrium 
configuration.  Assuming  that  =  f^{s^t)  is  a  dynamic  force  acting  on  the  thread  in  addition  to 

the  static  force  ,  namely  we  have  f  =  ,  then  Eq.  (10)  can  be  rewritten  as 


ds 


^  dv 

V  ds 


+  m- 


d^Y 


dt 


2  ’ 


(12) 


where  r  denotes  the  dynamic  deformation  from  the  static  equilibrium  position,  and  r  =  r  +  , 

where  is  the  solution  obtained  from  the  static  equilibrium  problem  Eq.  (11).  In  Eq.  (12),  we 

assumed  the  external  force  is  independent  of  the  tread  deformation.  However,  this  assumption  can 
be  relaxed  when  necessary  for  the  corresponding  applications.  Note  that  if  the  vibration  amplitude 
is  sufficiently  small,  then  Eq.  (12)  can  be  simplified  into  its  corresponding  linearized  form,  which 
will  be  discussed  in  the  next  section. 


4. 1. 1.2  Incremental  form 

Equation  (11)  results  in  a  two-point  nonlinear  boundary  value  problem.  To  solve  it,  an 

A  A  A 

incremental  form  of  the  thread  model  is  obtained  by  assuming  r  =  +  Ar,  T  =  7"^  + AT  ,  and 

+  Af ,  where  Fq  is  an  approximate  solution  of  Equation  (11)  obtained  from  a  prior  iteration 

A 

with  associated  tension  measure  and  loading  .  Using  an  important  relationship: 


-  (ir  /lEA  (  dlSY 

AT —  =  r— TT^  - 

ds  \ds 


the  incremental  form  of  Eq.  (11)  can  be  obtained  as 


d 

ds 


ds  y 


-GAr  =  e 


(13) 


(14) 
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where  B  is  called  (linearized)  stiffness  matrix  of  the  thread, 


^  FA  ^ 

l+frTX 

V  ^0  J 


(15) 


G  is  the  tangent  matrix  of  the  external  force, 


dr 


(16) 


and  e  is  the  residual  vector  obtained  at  the  last  iteration  step. 


^  f  0  ~T~ 

ds 


f 


dVr 


T  """0 
ds  j 


(17) 


Stiffness  matrix  B  in  Eq.  (16)  is  a  function  of  the  previous  tension,  Tq,  where  denotes  the 
tension  vector  T  obtained  at  the  previous  iteration  step,  and  is  the  generalized  tension 

-  EA  •  r 

obtained  at  the  previous  iteration  step,  ^ ,  where  =  Tq  .  Tangent  matrix  G  in  Eq. 

EA  +  T^ 

(17)  is  empty  when  f  is  independent  of  the  position  vector  r,  (for  example,  if  only  self  weight  of 
the  thread  is  considered),  and  it  becomes  non-empty  when  f  is  a  function  of  the  position  vector  r, 
(for  example,  if  the  thread  is  in  contact  with  other  deformable  objects.) 

Equation  (15)  is  linear  in  terms  of  Ar  .  Solving  Eq.  (15)  for  Ar  ,  an  updated  solution  for  r 
can  be  obtained  as 


r  =  Fq  +  Ar  , 


(18) 


Using  the  updated  r  ,  the  updated  tension  vector  T  can  be  obtained  as 


-  dr 
T  =  T 

ds  ’ 


(19) 


where 


^  ea-t 

T  = - and  T  =  EA 

EA  +  T 


dr 

\ds  J 


dr 

\ds  J 


-1 


(20) 
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Substituting  Eqs.  (20)  and  (21)  into  Eq.  (16)  the  stiffness  matrix  B  can  be  then  updated,  which 
can  be  then  used  for  next  iteration. 


Incremental  form  of  the  vibration  equation  Eq.  (12)  can  be  obtained  by  extending  the  above 
formulation  with  added  inertial  term.  As  a  result,  we  have 


—  B -  +m — ^  =  f^, 


ds 


(21) 


where  B  =  B(To) ,  =  T(rJ  ,  and  Tq  and  satisfy  the  static  equilibrium  equation: 


A. 

ds 


(  ^ 
T 

V 


ds 


+  f 


=  0. 


(22) 


It  is  assumed  here  that  vibration  amplitude  is  sufficiently  small,  so  the  stiffness  matrix  B  needs 
not  be  updated  during  the  vibration  analysis.  Therefore,  Eq.  (21)  becomes  a  linear  equation  for 
solving  Ar  .  Note  that  the  stiffness  matrix  B  in  Eqs.  (21)  is,  however,  not  constant  with  respect 
to  the  arc-length  coordinate,  and  in  general  we  have  B  =  B(^)  ,  which  is  determined  by  solving 
the  static  equilibrium  equation  (22). 


4. 1.1. 3  Weak  form  and  finite  element  discretization 

The  weak  form  of  Equation  (14)  is 


dSw^  _  dAr 


“'o  ds 


B 


ds 


-ds-^^  Svf^GArds  =  £  Svf^i^ds  -  £  f 


dSw^  dr. 


ds  ds 


-ds 


(23) 


where  Sw  denotes  virtual  displacement,  and  L  is  the  total  length  of  the  thread. 


A  finite  element  method  can  be  then  employed  to  discretize  Eq.  (23).  In  general,  assuming 
that  the  thread  shown  in  Fig.  1  is  divided  into  finite  elements,  each  element  has  two  nodes  at 

the  left  and  right  ends  with  coordinates,  s^  and  s^^^ ,  in  the  arc-length  coordinate  system  (where 
e  =  Also  assuming  that  0  =  0(^)  (^  e  [^^,^^+1])  is  the  selected  shape  function 

matrix  for  the  e-th  element,  then  for  element  e  we  have 


Ar  =  , 


(24) 
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where,  =|ArJ,Ar^^j|  is  the  vector  of  nodal  variable  perturbations  for  the  e-th  element. 
Then  the  equilibrium  equation  for  the  e-th  element  becomes 

(25) 

where 


f  = 

e  I  e 


d  d^ 


ds 


where 


ds 


(26) 

(27) 


Upon  assemble,  the  incremental  equation  for  the  residual  becomes 

KU  =  F,  (28) 

{T  T  T 

Ar^  ,  ,  •  •  • ,  Ar^  |  is  the  global  perturbation  vector  for  the  nodal  variables,  K  is 

the  global  stiffness  matrix,  and  F  is  the  global  residual  vector, 

K=A(k.)!uKiF  =  A(f,).  (29) 

Solving  Eq.  (28)  for  Ar^  ( e  =  1,2,*  •  + 1 ),  the  nodal  coordinate  vector  can  be  updated  as 


(30) 

where  is  the  nodal  coordinate  vector  of  element  e  obtained  in  the  last  step  and 
(e  =  l,2,“*,w^^  +1).  Using  the  updated  r^,  the  updated  tension  vector  T  for  element  e  can  be 
obtained  as 

T=fQ>u,  (31) 

e  e  e  e  ^  v  7 

where  ( e  =  1, 2,  •  •  • ,  + 1 ),  =  {rj  , 
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-  EA-T  (  i—f — 7:p — ^ —  \ 

’  and  7;  =  EA^u^O^  -l) • 


(32) 


Using  Eqs.  (31)  and  (32),  the  stiffness  matrix  for  element  e  can  be  updated,  as 


B  =7: 

e  e 


f  fA  ^ 

I  +  ^TTj 

rTi3  s  e 

\  J 


which  can  be  then  used  for  next  iteration  step. 


The  finite  element  equation  for  the  vibration  problem  ean  be  obtained  as 


KU  +  MU  =  F^ 


where  M  is  the  global  mass  matrix, 


M  =  A(m,) 


and 


(33) 


(34) 


(35) 


(36) 


4.1.2  Numerical  implementation 

4. 1.2. 1  Finite  element  discretization  with  linear  shape  functions 

With  linear  interpolation,  the  ineremental  elementary  nodal  displaeement  veetor,  ,  is 
assumed  to  have  the  following  form: 


(37) 


where  Ar.  =  {ax.,aj7,.,azJ  is  the  incremental  nodal  displacement  vector  for  the  i-th  node.  The 
shape  function  matrix,  ,  is  then  assumed  to  have  the  following  form; 


0  0  ^2  0  0 
0  0  0  ^2  0 
0  0  0  0  ^2 


(38) 
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where, 


(I)^=Y{sm-s),  <2^2  =-^(^-‘y/).and  J,.  -5,.  for  ^  e 


(39) 


4. 1.2.2  Finite  element  discretization  with  second-order  polynomials 

With  the  seeond-order  polynomial  interpolation,  nodal  displaeement,  is  assumed  to  have 
the  following  form: 


u 


e 


(40) 


where  a  mid-node  with  the  incremental  displacement,  added 

in  the  second-order  element.  The  shape  function  matrix,  O  ,  is  then  assumed  to  have  the 
following  form: 


3>„  = 


>1  0  0 
0  0 
0  0 


^2  0  0  ^3  0  0 

0  ^2  0  0  ^3  0 

0  0  ^2  0  0  ^3 


where, 


(41) 


,  ,  3s  2s^  ,  As  4^^ 

®=1 +  — T’  f2= - T 

‘  L  r  ^  L  r 


,  S  2s^  „  r/^  n 

<2)3  = -— -h ,  for  5  e  [0,Z] 


(42) 


With  the  assumption  that  the  external  foree  is  independent  of  fiber  position  (namely,  G  =  0), 
the  element  stiffness  matrix  ean  be  expressed  as 


=  f  O^BOds  =  r 

Jo  Jo 


O'^B^  ds  +  1 , 0  ‘  B^  ds 


'T-, 


(43) 


For  further  numerical  implementation,  The  matrix  O  in  Eq.  (28)  is  expanded  below: 
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^1^11^1  ^1^12^1  ^1^13^1  ^1^11^2  ^1^12^2  ^1^13^2  ^1^11^3  ^1^12^3  ^1^13^3 

^1^22^1  ^1^23^1  ^1^21^2  ^1^22^2  ^1^23^2  ^1^21^3  ^1^22^3  ^1^23^3 

^1^33^1  ^1^31^2  ^1^32^2  ^1^33^2  ^1^31^3  ^1^32^3  ^1^33^3 

^2^11^2  ^2^12^2  ^2^13^2  ^2^11^3  ^2^12^3  ^2^13^3 

^2^22^2  ^2^23^2  ^2^21^3  ^2^22^3  ^2^23^3 

^2^33^2  ^2^31^3  ^2^32^3  ^2^33^3 

^3^11^3  ^3^12^3  ^3^13^3 

^3^22^3  ^3^23^3 

^3^33^3 

In  the  element  level,  the  residual  veetor  is  expressed  as  follows: 

fo'^r .  (44) 

The  seeond  term  in  Eq.  (44)  ean  be  further  expressed  as  follows 

r  ^or,’ 

jiri'-rol-.  .  (45) 

[  -?;>•,' 

where  and  r^  are  determined  in  domain  ^g[0,Z/2]  ,  and  and  are  determined  in 
domain  ^  g  [Z/2,Z] . 

4.1 .3  Extended  formula  for  multi-thread  systems 

To  form  the  initial  eonfiguration  of  a  web,  multiple  threads  are  plaeed  in  the  Cartesian 
eoordinate  system  (O-XYZ  in  Fig.  5-a)  in  whieh  the  eonneeted  nodes  among  threads  share  the 
same  eoordinate.  At  the  initial  state,  the  eonfiguration  of  the  web  is  estimated  and  a  seareh  for 
the  equilibrium  is  followed. 

In  the  fabrie  web,  eaeh  thread  is  diseretized  in  its  own  are-length  eoordinate,  5.  The 
ineremental  equation  for  eaeh  thread  ean  be  obtained  using  the  same  formulations  as  Eq.  (28)  or 
(34).  The  ineremental  equations  for  individual  threads  are  then  assembled  together  to  form  the 
global  stiffness  matrix,  K,  and  the  global  residual  veetor,  F,  for  the  web  system.  The  assembly 
proeedure  is  same  as  the  stiffness  and  foree  veetor  assembly  proeedure  in  any  eommereial  FEM 
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code.  With  the  assembled  finite  element  equation  for  the  whole  web,  the  incremental 
displacement,  U,  ean  be  obtained  to  update  the  configuration  of  the  web.  By  this  means,  the 
statie  equilibrium  of  the  web  system  under  any  external  load  ean  be  achieved.  The  solution 
proeedure  for  the  vibration  problem  is  similar  to  the  static  problem  and  is  not  repeated  here. 

4.1.4  Examples 

An  analysis  eode  has  been  then  developed  to  simulate  the  large  nonlinear  deformation  and  high- 
fidelity  stress  distribution  for  the  fabric  system  with  arbitrary  fiber  configurations  and  under  statie 
and  dynamic  loads.  This  code  is  being  extended  to  deal  with  more  general  textile  fabrics  with 
non-linear  geometry  effeets  and  to  be  embedded  in  the  meso-structural  model  for  the  textile 
eomposites. 

Figure  6  demonstrates  the  new  capability  developed  with  solving  a  statie  equilibrium 
problem  of  a  three-dimensional  spider  net-like  fiber  net.  Figure  6-a  illustrates  the  geometry  of  the 
spider  net,  while  Fig.  6-b  shows  the  resultant  static  deformation  predieted  using  the  formulation 
developed.  Note  that  carbon  fiber  braided  spider  net  may  be  used  as  the  reinforcement  in 
composite  structures.  The  analysis  code  has  been  developed  to  simulate  the  deformation  and 
stress  for  any  earbon  fiber  net  under  given  static  and  dynamic  loads. 


Figure  6.  The  shape  of  a  spider  net  and  the  corresponding  deformation. 

Figure  7  shows  another  example  that  an  out-of-plane  vertical  load,  F  =  50,884  N,  is  applied 
to  the  spider  net  in  the  negative  z-direction.  The  deformed  net  and  the  tension  force  in  the  net  are 
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illustrated  in  Fig.  7-a  and  7-b,  which  shows  a  maximum  tension,  37,225  N,  in  the  net  under  the 
given  maximum  impact  load.  The  red  color  in  Fig.  7-b  indicates  the  maximum  tension,  while  the 
other  colors  indicate  other  tension  levels  as  shown  in  the  color  map. 
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(a)  Deformation  shape 


(b)  Tension  distribution 


Figure  7.  A  concentrated  force,  F=50,884  N,  is  applied  at  the  red  point  in  negative  z- 
direction.  a)  Spider  net  deformation  under  the  force  F.  b)  The  maximum  tension  force  in 
the  net  is  shown  in  red  color  with  an  amount  of  37,225  N. 


Figures  8  and  9  illustrate  an  example  for  demonstrating  the  new  fiber  modeling  capability 
developed  to  simulate  a  tennis  racket  during  tennis  ball-racket  interaction.  The  objective  of  the 
simulation  example  is  to  find  out  how  a  pretension  can  be  determined  for  certain  design  purposes 
related  to  racket  performance  with  a  given  tennis  racket  configuration.  In  this  example,  a  tennis 
racket  shown  in  Fig.  8  is  considered,  which  has  16x19  strings  with  a  diameter  of  1.24  mm.  The 
tensile  modulus  of  the  string  is  £  =  90  GPa  with  a  recommended  pretension  between  55  N  and  62 
N.  Numerical  tests  have  been  done  for  various  pretensions  and  forced  out-plane  displacements 
due  to  tennis  ball-racket  interactions.  Figure  9-a  shows  an  example  result  for  predicted  tension 
distribution  in  the  racket  strings  under  a  given  pretension  and  forced  displacement,  while  Fig.  9-b 
illustrates  the  out-plane  stiffness-displacement  relationships  under  different  pretensions.  This  is 
seen  that  out-plane  stiffness  of  the  tennis  racket  is  non-linear  with  respect  to  the  out-plane 
displacement  and  will  increase  nonlinearly  with  the  pretension.  This  kind  of  results,  with  other 
available  results  from  the  simulation,  for  example,  failure  of  the  strings,  could  be  very  useful  for 
designing  an  advance  tennis  racket.  More  discussions  regarding  how  to  use  the  new  capability  for 
real  tennis  racket  design  problem  is  however  omitted  in  this  report. 
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Figure  8.  A  tennis  racket  simulation  problem 


Stiffness  vs.  Displacement  (same  pretension  fer  a  line) 


a)  Example  tension  distribution  b)  Predieted  raeket  out-plane  stiffness  vs. 

maximum  displaeement  and  pretension 

Figure  9.  Example  results  using  the  new  fiber  model 


The  above  examples  have  demonstrated  the  effeetiveness  of  the  new  fiber  model  and 
assoeiated  simulation  eapabilities  developed  for  designing  a  fiber  or  eable  net  system  with  high 
fidelity  under  arbitrary  large  deformation. 

4.2  High  Fidelity  Meso-Structural  Modeling  of  Textile  Composite 

In  this  task,  we  will  develop  a  high-fidelity  meso-struetural  model  for  general  textile  eomposite, 
whieh  will  inelude  fiber/fabries  and  filling/matrix  material.  This  meso-struetural  model  will  allow 
us  to  eonduet  a  detailed  loealized  analysis  for  the  meso-behaviors  of  the  eomposite,  whieh 
inelude  nonlinear  deformation  and  strength  of  individual  fibers  in  the  eomposite. 
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4.2.1  Modeling  of  fiber-fiber  interaction 

In  the  previous  eable  net  models  (e.g.,  in  Fig.  6),  it  is  assumed  the  eross  threads  share  the  same 
node  at  the  erossing  point,  whieh  implies:  1)  eables  are  tied  together  with  a  knot  and  relative 
movement  (e.g.,  slip)  of  erossing  threads  is  not  allowed  at  the  knot,  2)  effeet  of  thread  thiekness 
(diameter)  is  ignored,  and  3)  most  importantly,  eontaet  foree  between  the  threads  ean  not  be 
aeeurately  predieted.  In  the  eurrent  researeh,  the  nonlinear  fiber  model  has  been  extended  to 
eonsider  thread  to  thread  eontaet  and  effeet  of  eross-seetion  geometry  of  the  threads.  Figure  10 
illustrates  the  new  fiber- fiber  interaetion  model  (Fig.  10-b)  with  a  eomparison  with  the  previous 
eable  net  model  (Fig.  10-a)  at  the  node. 


a)  Previous  eable  net  model  b)  New  fiber-fiber  interaetion  model 


Figure  10.  Previous  and  new  fiber-fiber  interactions  models 

First,  eaeh  thread  is  eonsidered  deformable  along  the  radius  direetion  with  a  diameter  and 
stiffness  eoeffieient.  The  aetion-reaetion  foree  between  the  threads  is  then  determined  by  the 
relative  deformation  of  the  eontaeting  treads.  The  deformation  of  eaeh  thread  is  also  a  funetion  of 
the  eontaet  foree.  An  iterative  proeess  has  therefore  been  developed  to  obtain  the  eonverged 
results  for  both  the  eontaet  forees  and  threads  deformations. 

Figure  11  and  Table  4  eompare  the  simulation  results  obtained  with  the  two  different  fiber- 
fiber  interaetion  models  deseribed  in  Fig.  10.  Figure  1 1  illustrates  the  deformations  (amplified  by 
three  times)  and  tension  distributions  in  the  woven  fabries,  while  Table  4  shows  the  out-plane 
stiffness  and  maximum  and  minimum  tension  values.  Exaetly  the  same  boundary  eonditions 
were  used  in  the  two  simulations,  whieh  fix  the  two  ends  of  eaeh  fiber.  A  eoneentrated  foree, 
F=5001bs,  is  applied  at  the  eenter  of  eaeh  woven  fabrie.  As  shown  in  Fig.  1 1  and  Table  4,  while 
higher  tension  is  predieted  with  the  new  fabrie  model,  the  new  model  tends  to  distribute  the 


23 


tension  force  more  evenly  in  the  fibers.  As  shown  in  Table  4,  the  new  fabric  model  also  results  in 
a  higher  out-plane  stiffness  under  the  given  force.  Figure  11 -a  also  shows  that  the  maximum 
tension  in  the  middle  weft  if  the  previous  fabric  model  is  used,  while  Fig.  11-b  shows  that  the 
middle  weft  will  have  lower  tension  compared  with  the  other  two  wefts.  Physical  explanation  for 
this  is  that  the  tension  in  the  middle  weft  has  been  partially  relieved  due  to  the  applied  force  on 
the  middle  weft.  This  phenomenon  can  not  be  predicted  with  the  previous  fabric  model.  Therefore, 
the  new  fiber-fiber  contact  model  is  more  accurate  then  the  previous  model. 


a)  Previous  cable  net  model 


b)  New  fiber-fiber  interaction  model 


Figure  11.  Comparison  of  previous  and  new  fiber-fiber  interaction  models 


Table  4.  Comparison  of  predicted  out-plane  stiffness  and  tension 


Out-plane 

Stif&ess 

Maximum 

Tension 

Minimum 

Tension 

Tension 

Difference 

Previous  Model 

1,020  lb/ in 

3,470  lb 

2,060  lb 

1,4101b 

New  Model 

1,429  Ib/in 

4,230  lb 

3,550  lb 

6801b 
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Figure  12  illustrates  that  the  new  fiber-fiber  eontaet  model  ean  be  further  extended  to  deal 
with  real  woven  and  braided  fabries  design  problems  with  the  design  variables,  sueh  as:  yam 
width  (Fig.  12-b  and  e),  woven  angle  (Fig.  12-d),  relative  size  (Fig.  12-e),  and  woven  forees  (Fig. 
12-f).  The  preliminary  results  obtained  in  this  report  have  demonstrated  that  the  eapability 
developed  ean  be  extended  to  evaluate  widely  different  designs  of  woven  eomposites  for  the 
purpose  of  improving  stmetural  performanee  and  redueing  weight.  By  applying  an  optimization 
proeess  as  diseussed  in  Seetion  2.1.1,  the  teehnique  ean  be  further  extended  to  find  the  optimum 
arrangement  of  the  woven  fibers  in  the  eomposite  for  improved  performanee  subjeet  to  the  given 
weight,  loading  eonditions,  manufaeturing  eost,  and  other  eonstraints. 


c)  Change  Yarn  Width-2 


a)  Nominal  Design 


b)  Change  Yarn  Width- 1 


d)  Change  Woven  Angle 


e)  Change  Relative  Sizes 


f)  Change  Woven  Forces 


Figure  12.  Example  designs  of  a  woven  composite  material  (refer  to  [16]). 


4.2.2  Development  of  an  innovative  composite  concept 

We  are  utilizing  the  new  Funetion-Oriented  Material  Design  (FOMD)  methodologies  to  explore 
innovative  stmetural  and  material  designs.  One  example  new  eoneept  material  developed  is 
termed  Biomimetic  Tendon-Reinforced  (BTR)  eomposite  material.  Figure  13  deseribes  a  design 
proeess,  whieh  demonstrates  the  feasibility  of  FOMD  to  lay  out  innovative  eoneepts  of  eomposite 
materials.  The  example  problem  is  to  design  a  multi-layer  eomposite  material  that  features 
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functional  outer  and  inner  covers  and  a  typical  thickness  as  shown  in  Fig.  13-a.  The  middle  layer 
between  the  outer  and  inner  covers  is  considered  to  be  a  composite  of  two  different  materials:  one 
is  relatively  stiff  and  forms  the  core  structure  (e.g.,  fiber  reinforcements)  for  the  major 
reinforcement;  the  other  material  is  soft  and  aets  as  the  matrix  for  the  composite  (which  could  be, 
for  example,  a  polymer).  Both  materials  are  considered  as  fiinetional  materials.  For  instance,  the 
fiber  material  may  carry  the  major  loads  applied  to  the  structure,  while  the  matrix  material  could 
be  an  insulative  material  for  the  purpose  of  thermal,  vibration,  and/or  acoustie  managements.  The 
objective  considered  in  this  example  problem  is  to  maximize  the  out-of-plane  stiffness  and 
strength  due  to  an  out-of-plane  load  by  minimizing  the  total  strain  energy  stored  in  the  strueture. 
The  constraint  function  is  the  total  amount  of  material  used  to  build  the  strueture  between  the 
outer  and  inner  covers.  The  design  problem  is  defined  as  shown  in  Fig.  13-a,  and  the  optimal 
design  for  the  design  problem  is  obtained  using  the  FOMD  proeess,  which  is  shown  in  Fig.  13-b. 
Figure  14  illustrates  that  the  two-dimensional  design  in  Fig.  13-b  ean  be  extended  to  a  three- 
dimensional  configuration. 


Optional  inner  cover 

Joint 


Cover 


Design  Domain 

a )  Design  P  ro  b  le  m 

S  tuffer 


O  ptional  outer  cover 


Functional  materials 

(carbon  fiber  and  epoxy  matrix) 


Fiber 


M  atrix 


b)  Optimum  structure  obtained 


Figure  13.  Design  problem  for  back  plate. 
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d)  Matrix 


e)  Covers 


Figure  14.  From  two-dimensional  design  to  a  three-dimensional  configuration. 

The  new  eoneept  material  is  termed  Biomimetic  Tendon-Reinforced  (BTR)  eomposite 
material.  As  shown  in  Fig.  15,  the  BTR  eomposite  ineludes  five  fundamental  eomponents: 
tendons/museles  (represented  by  fiber  eables  and/or  tendon  aetuators),  ribs/bones  (represented  by 
metallie  or  eeramie  staffers  and  struts),  joints  (ineluding  knots),  flesh  (represented  by  filling 
polymers,  foams,  ete.),  and  skins  (represented  by  woven  eomposite  layers  or  other  thin  eovering 
materials).  Various  BTR  material  eoneepts  have  been  developed  at  U  of  M  in  eollaborating  with 
MKP  Ine.  for  potentially  different  applieations. 


27 


Biomimetic  Tendon-Reinforced  (BTR)  Material 
(©  MKP  Patent  Pending  ) 

Five  fundamental  components: 

•  Tendons/muscles  -  fiber 
cables  and/or  actuators 

•  Ribs/bones  -  metallic  or 
ceramic  stuffers  and  struts 

•  Joints  -  mechanical  joints, 
fasteners,  knots,  etc. 

•  Skins  -  woven  composite 
layers  or  other  covering 
materials 

•  Flesh  -  filling  polymers,  foams, 
etc. 

Figure  15.  BTR  composite  material  concept  developed  at  MKP  Inc. 


Figure  16  illustrates  two  example  BTR  eoneepts  and  assoeiated  fabrieation  proeesses. 
Figure  16-a  shows  a  BTR  material  with  bead  stuffers,  earbon  fiber/epoxy  skin,  and  Kevlar 
tendons,  while  Fig.  16-b  shows  one  with  rod  stuffers,  Kevlar  fiber/epoxy  skin,  and  Kevlar 
tendons.  These  materials  are  ultra-lightweight  and  ean  be  fabrieated  using  the  weaving  proeesses 
shown  in  Figs.  16-a  and  b.  The  BTR  eoneept  will  result  in  a  new  and  innovative  kind  of 
eomposite  materials:  unlike  traditional  woven  materials,  in  whieh  the  fibers  are  nearly  evenly 
distributed  in  one  plane  in  the  matrix  material,  the  new  material  will  be  reinforeed  by  alloeating 
eoneentrated  fibers,  sueh  as  fiber  ropes,  along  load  paths  so  as  to  inerease  transverse  stiffness. 
MKP  has  filed  a  patent  for  this  new  BTR  material  eoneept. 


a)  BTR  with  bead  stuffer  b)  BTR  with  rod  stuffer 

Figure  16.  Example  BTR  materials  and  fabrication  processes  developed  at  MKP  Inc. 
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We  are  exploring  the  use  of  the  above  eomposite  material  eoneept  for  various  applieations. 
Several  example  BTR  prototypes  have  been  developed  at  MKP,  as  shown  in  Fig.  17.  Figure  17 
illustrates  three  BTR  prototypes  (Fig.  17  b-d)  previously  developed  at  MKP.  We  will  seleet  one 
eonfiguration  of  the  eoneept  designs,  through  virtual  prototyping  and  experiments,  for  a  baek 
plate  in  a  eomposite  armor.  We  will  further  design  and  fabrieate  a  mold  with  the  fixtures  for 
fabrieating  the  baek  plate  in  a  eost-effeetive  and  effieient  way.  With  the  new  mold,  a  number  of 
sample  baek  plates  will  be  made,  whieh  will  be  used  for  basie  meehanieal  tests  and  possible 
ballistie  tests. 


Figure  17.  Options  for  composite  plate  design:  a)  Composite  plate  with  multi-layered  rod 
stuffers;  b)  Composite  plate  with  rod  stuffers;  c)  Composite  plate  with  bead  stuffers:  d) 

Conventional  laminate  composite 


The  potential  benefits  of  using  the  BTR  eomposite  have  been  investigated.  A  two- 
dimensional  eoneept  BTR  eomposite  is  shown  in  Fig.  18-a,  whieh  is  eomposed  of  glass  fiber 
tendons  and  alumina  stuffers.  For  eomparison,  a  eonventional  S2  glass  laminate  is  also  modeled 
with  the  same  size,  whieh  is  shown  in  Fig.  18-b.  The  elastie  property  and  strength  of  the  S2 
fiber/epoxy  material  and  alumina  material  are  listed  in  Table  5. 


(a)  BTR  strueture 


Epoxy  matrix 


S2  glass  fiber 


(b)  Glass  fiber  laminate  strueture 


Figure  18.  Finite  element  models  for  BTR  and  traditional  laminate  composites. 


29 


Table  5.  Material  properties  used  for  composite  analysis 


Glass  fiber 
with  epoxy 

31astie  Properties 

£,(GPa) 

(GPa) 

G,(GPa) 

S2/3501-6 

55.9 

18 

0.27 

0.27 

7.6 

Glass  fiber 
with  epoxy 

Composite  Strengt 

1 

5/  (GPa) 

(GPa) 

(GPa) 

5/  (GPa) 

*^12  (GPa) 

523(GPa) 

S2/3501-6 

1.79 

0.96 

0.06 

0.16 

0.11 

0.083 

Glass  fiber 
with  epoxy 

Fraeture  Properties 

Density 

Gi,  (JW) 

P  (g/cc) 

S2/3501-6 

160 

1.99 

stuffer 

Elastie  Properties 

Compressive  Strength 

Density 

E  (GPa) 

V 

(GPa) 

P  (g/cc) 

Alumina 

(99%) 

370 

0.22 

3.0 

3.96 

{E^  is  fiber  axial  modulus,  is  fiber  transverse  modulus,  is  axial  Poisson’s  ratio,  is 
transverse  Poisson’s  ratio,  and  is  axial  shear  modulus.  iS/  is  fiber  axial  tensile  strength,  is 
fiber/epoxy  axial  eompression  strength,  S2  is  fiber/epoxy  transverse  tensile  strength,  82^  is 
fiber/epoxy  transverse  eompression  strength,  S^2  fiber/epoxy  in-plane  shear  strength,  and  5^23 
is  fiber/epoxy  out-of-plane  shear  strength.) 

In  the  eurrent  simulation,  a  blast  load  is  represented  as  an  overpressure  on  the  eomposite. 
The  two-dimensional  eomposite  panel  in  Fig.  19  has  a  dimension  36  em  -  3  em  in  length  and 
thiekness.  The  panel  is  elamped  at  both  edges  along  the  length.  A  peak  overpressure,  112  MPa, 
is  reaehed  at  the  top  eenter  of  the  panel  with  an  8-kg  TNT  explosive  eharged  50  em  above  the  top 
eenter  of  the  panel.  The  pressure  was  assumed  to  deeay  exponentially  following  the  peak  load. 
The  struetural  responses  of  the  eomposite  panels  with  the  peak  overpressure  are  analyzed  and 
shown  in  Fig.  19. 


(b)  BTR  eomposite 

Figure  19.  Comparison  of  BTR  composite  with  traditional  laminate  composite  under  a 

blast  load  (tensile  stress  in  fibers.) 
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The  appropriate  failure  modes  of  the  two  eomposites  are  summarized  in  Table  6.  The 
weight  and  meehanieal  properties  of  the  BTR  eomposite  panel  and  laminate  eomposite  panel  are 
eompared  in  Table  7.  It  is  seen  that  the  BTR  eomposite  has  superior  strength  (approximately  2.5 
times)  eompared  to  the  eonventional  laminate  panel  under  the  blast  loading.  The  BTR  panel  is 
also  mueh  lighter  than  the  laminate  panel.  However,  the  out-of-plane  rigidity  of  the  BTR  panel  is 
redueed  to  55%  of  the  laminate  eomposite  panel. 


Table  6.  Failure  modes  for  BTR  and  laminate  composite  panels  under  blast 


Failure  mode 

BTR  composite 

Laminate  composite 

Fiber  failure 

Tensile  stress  in  fibers 

Tensile/compressive  stress  in  fibers 

Matrix/stuffer 

failure 

Staffer  compressive 
stress 

Transverse  normal  stresses  and  shear  stresses 

Delamination 

Not  applicable 

Through  the  thickness  normal  stress  and  applicable 
shear  stresses 

Table  7.  The  comparison  of  BTR  composite  and  laminate  composite 


Weight 

(normalized) 

Out-of-plane  rigidity 
(normalized) 

Strength* 

(normalized) 

BTR 

68% 

55% 

245  % 

Laminate 

100  % 

100  % 

100  % 

(*  The  BTR  strength  index  is  ealeulated  as  the  maximum  fiber  tensile  stress  divided  by  fiber 
axial  tensile  strength;  the  laminate  strength  index  is  ealeulated  using  the  equation  in  Yen  and 
Jones  [17]. ) 


4.2.3  Modeling  of  fiber  interactions  with  other  materials  in  a  composite 

The  true  nonlinear  fiber  model  has  been  extended  to  eonsider  interaetion  with  other  materials, 
for  example,  stuffers  and  matrix,  in  a  eomposite.  The  resultant  eapability  has  been  utilized  to 
simulate  response  of  a  BTR  eomposite  under  a  blast  loading.  Figure  20-a  shows  the  predieted 
displaeement  and  tensile  stress  in  the  BTR  strueture  under  the  same  blast  loading  as  applied  in 
Fig.  20-b,  while  Fig.  20-b  shows  the  results  obtained  using  a  eonventional  finite  element  model. 
The  predieted  maximum  tensile  stress  in  the  fiber  tendon  is  385  MPa  using  the  true  nonlinear 
model,  whieh  is  larger  than  the  maximum  axial  tensile  stress  (307  MPa)  obtained  from  the 
eonventional  FE  model.  Sinee  the  non-eompressibility  of  the  fiber  tendon  is  eonsidered  in  the 
eurrent  model,  the  result  in  Fig.  20-a  is  supposed  to  be  more  aeeurate  than  the  FEA  result.  The 
tensile  stress  distribution  obtained  from  the  true  nonlinear  eable  model  is  also  mueh  smoother 
than  the  result  from  eonventional  FE  model,  whieh  indieates  that  the  nonlinear  eable  model  ean 
eapture  the  instant  stress  wave  propagation  in  the  eable  network  under  blast.  The  true  nonlinear 
eable  model  will  be  further  extended  and  eorrelate  with  the  experimental  results  and  then 
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implemented  into  optimal  design  system  for  designing  advaneed  BTR  stmeture  with  mueh 
improved  performanee. 


385.17 

330.14 

275.12 

220.09 

165.07 

110.05 

55.024 

0 


(a)  Result  obtained  using  the  true  nonlinear  eable  model 


>2.19e+02 
<219e+02 
<131e+Q2 
<  439e+01 
<-4.36e+01 
<-1.31e+02 
<-2.ige+02 
-:-3.06e+02 


Max  =  3.07e+02 
Min  =-3.06e+02 


(b)  Result  obtained  using  the  eonventional  finite  element  model 


Figure  20.  Comparison  of  tensile  stress  (MPa)  in  fiber  tendons  under  blast  loading  using 
different  simulation  models. 


The  BTR  eomposite  eoneept  features  ultra-light  weight  and  outstanding  out-of-plane 
stiffness/strength,  whieh  ean  also  provide  superior  proteetion  against  shoek  wave  in  a  blast  event. 
Note  that  aeeording  to  our  knowledge,  baek  plate  will  have  signifieant  effeet  on  armor 
performanee.  Our  baek  plate  with  BTR  eomposite  eoneept  will  be  designed  to  have  maximum 
bending  stiffness  and  strength  with  a  minimum  weight  for  given  armor  applieations.  The 
preliminary  results  shown  in  this  report  ean  be  further  extended  to  deal  with  other  problems  sueh 
as  eompression  strength,  buekling  stability,  and  impaet  resistanee,  whieh  will  be  major  tasks  in 
next  phase  of  this  researeh  program.  Moreover,  multiple  objeetive  funetions  ean  be  eonsidered 
for  designing  a  eomposite  stmeture  with  multiple  roles  in  an  advaneed  vehiele  system. 

The  non-linear  fiber  model  is  being  further  extended  to  eonsider  fiber-matrix  interaetion  in  a 
general  eomposite  material.  Figure  21  illustrates  sueh  geometry  model  has  been  developed  with  a 
foeus  on  predieting  meso-behavior  of  a  eomposite  and  its  delamination  and  de -bonding 
phenomena  as  well  as  assoeiated  damage  meehanism. 
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Figure  21.  An  example  model  of  fiber-matrix  interaction. 

In  summary,  we  have  extended  our  nonlinear  fiber  model  to  general  three-dimensional  fabrie 
system  with  various  eonfigurations  for  simulating  praetieal  textile  fabries  and  eomposite 
materials.  The  formulation  will  aeeount  for  geometry  nonlinearity  and  arbitrary  large  deformation 
of  the  textile  fabries.  This  results  in  a  unique  eapability  for  simulating  textile  eomposite  materials, 
espeeially  for  the  damage  assessment.  The  new  eapability  will  be  further  extended  and  integrated 
with  the  other  advaneed  eapabilities  being  developed  for  designing  truly  optimum  eomposite  for 
various  applieations. 

4.2.4  Example  results 

4. 2. 4. 1  Simulating  a  simple  multi-fiber  interaction  model 

Figure  22  illustrates  a  simple  model  for  a  fiber- fiber  eontaet  problem  with  two  erossing-over 
threads.  As  shown  in  Fig.  22,  eaeh  thread  is  fixed  at  the  both  ends  with  a  given  pretension,  whieh 
is  defined  by  an  initial  tensile  strain,  (where  i  =  1,2).  The  thread-to-thread  interaetion  results 

in  a  eontaet  foree  between  the  two  threads,  whieh  is  a  funetion  of  pre-seleeted  design  parameters, 
sueh  as  thread  diameter  D  (for  the  both  threads)  and  pre-tensions  (defined  by  initial  tensile  strains 
^^and  ^2  ill  Fig-  22).  Three  examples  have  been  simulated  with:  1)  D=2.0,  e^=e2=0.94% ,  2) 

D=2.0,  e^  =0.94%,  ^2  =2.0%  ,  and  3)  Z)=1.2,  =62  =0.94%  .  The  same  tensile  modulus, 

EA=2.5q5,  is  used  for  all  the  three  examples. 


Figure  22.  Simulation  model  of  a  simple  fiber-fiber  contact  problem 

Figure  23  shows  the  results  obtained  using  the  new  fiber-fiber  eontaet  model  for  the  three 
examples.  It  is  seen  that  inereasing  the  initial  strain  in  thread  #2  by  113%  results  in  a  23% 
inerease  of  the  eontaet  foree  (eomparing  Fig.  23-b  with  Fig.  23-a),  while  redueing  thread  diameter 
by  40%  will  reduee  the  eontaet  foree  by  81%.  With  this  kind  analyses  enabled  by  the  new 
eapability  developed,  we  ean  design  a  fabrie  system  with  neeessary  detailed  information  for 
various  design  purposes,  sueh  as  redueing  potential  fiber  failures  in  a  fabrie  system. 


a)  D=2.0,  61=62=0.94% 


b)  D=2.0,  61=0.94%  62=2% 


c)  D=1.2,  61=62=0.94% 


Figure  23.  Example  results  for  using  different  design  parameters 
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4. 2. 4. 2  Simulating  a  braided  fabric 

The  second  example  is  to  demonstrate  the  new  capability  developed  can  be  used  to  simulate  a 
very  general  fabric  system.  A  simulation  model  with  a  warp  fiber  and  three  weft  fibers  is  shown 
in  Fig.  24.  As  shown  in  Fig.  24,  each  fiber  in  the  model  is  fixed  at  the  both  ends,  and  a  2.04% 
pre-tension  strain  is  applied  initially  to  all  the  fibers.  Figure  25  and  Table  5  illustrate  example 
simulation  results  obtained  using  the  simulation  code  for  four  different  cases:  a)  nominal  case, 
EA=2.5q5,  D=2.0;  b)  doubled  tensile  modulus,  namely,  EA=5.0e5,  D=2.0;  c)  reduced  fiber 
diameter,  EA=2.5e5,  D=1.0;  and  d)  only  the  warp  tensile  modulus  is  doubled,  namely,  EA=5.0e5 
for  the  warp,  EA=2.5e5  for  wefts,  and  £)=2.0.  Figure  25  shows  that  tension  distribution  in  the 
fabric  system  becomes  very  different  under  different  design  parameters.  For  example,  reducing 
fiber  diameters  of  the  fibers  results  in  tension  to  be  more  evenly  distributed  in  the  braided  fibers 
(comparing  Fig.  25-c  with  Fig.  25-a),  while  using  different  fibers  for  the  ward  and  wefts  (e.  g., 
doubling  the  warp  tensile  modulus  as  shown  in  Fig.  25 -d)  results  in  more  unevenly  distributed 
tension  in  the  fabric  system.  Table  5  further  compares  the  actual  tension  values  for  the  four  cases, 
which  include  maximum  and  minimum  tensions  and  difference  between  them.  It  is  seen  that 
doubling  the  tensile  modulus  for  all  the  fibers  will  almost  double  the  tension  in  the  braided  fibers, 
while  doubling  only  the  warp  tensile  modulus  will  results  in  a  23%  increase  in  the  maximum 
tension  and  31%  increase  in  the  minimum  tension.  Again,  these  kind  results  can  be  used  to 
design  a  braided  fabric  with  desired  material  properties  and  damage  tolerance. 


Figure  24.  Deformation  of  braided  fibers  under  a  2.04%  pretension  strain. 
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b)  EA=5.0e5  for  all 


Figure  25.  Comparison  of  predicted  tensions  for  different  parameters 
(EA:  tensile  modulus;  D:  diameter  of  the  fibers). 


Table  5.  Comparison  of  tension  values  (N) 


Maximum 

Tension 

Minimum 

Tension 

Tension 

Differenee 

a)  EA=2.5e5,  D=2.0 

4,429 

3,348 

1,081 

b)EA=5.0e5,D=2.0 

8,191 

6,029 

2,162 

c)EA=2.5e5,D=l.O 

3,408 

2,315 

1,093 

d)  EA=5.0e5  for  warp,  D=2.0 

5,456 

4,311 

1,145 

4.2.43  Simulation  of  a  waiting  material  concept 

Figure  26  illustrates  an  example  eomposite  strueture,  whieh  is  made  of  earbon  fibers  and 
aluminum  stuffers.  This  strueture  represents  a  novel  baekplate  eomposite  eoneept  for  an 
advaneed  armor,  whieh  is  also  used  to  illustrate  the  eoneept  of  multi-stage  stability  of  the  eoneept 
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structure.  The  key  idea  of  the  sample  composite  structure  is  the  embedded  woven  fiber  ropes. 
When  a  load  carrying  thread  is  broken,  the  neighboring  threads  can  act  as  the  safety  members  to 
reserve  the  integrity  of  the  whole  structure  provided  the  threads  are  properly  plaeed.  As  shown  in 
Fig.  26,  five  metallie  beads  are  utilized  as  the  stuffers  in  a  braiding  proeess  to  form  a  woven 
lattiee  eomposite.  The  integrity  of  the  eomposite  strueture  is  supported  by  the  pre-tensioned 
fabrie  ropes.  When  a  rigid  objeet  is  impaeted  on  the  eomposite,  the  deformation  of  the  strueture 
and  the  eorresponding  tension  foree  in  the  rope  ean  be  obtained  by  using  the  nonlinear  analysis 
model  developed  by  us. 


Figure  26.  The  sample  composite  structure  for  multi-stage  stability  illustration 

Figure  27  illustrates  the  basie  eoneept  of  the  multi-stage  stability  in  the  eomposite  strueture. 
The  maximum  permissible  tensile  foree  in  the  threads  is  3,800  N.  In  Fig.  27-a,  the  flying  objeet 
hits  the  eomposite  strueture,  the  maximum  defleetion  of  the  eomposite  strueture  beeomes  2.5  mm, 
while  the  magnitude  of  the  tensions  in  the  threads  is  shown  in  the  eolormap.  It  is  seen  that  the 
tension  in  the  master  thread  (shown  in  red)  is  elose  to  the  strength  limit,  and  the  neighboring 
thread  (shown  in  green)  is  going  to  take  effect  in  the  next  stability  stage.  In  Fig.  27-b,  the 
stability  stage  B  reaehes  its  limit,  the  red  thread  is  going  to  break,  while  the  eyan  neighboring 
thread  is  going  to  aet  in  stability  stage  C.  Figure  27-c  shows  the  stability  stage  C.  It  is  seen  that 
the  central  metal  staffer  is  separated  from  the  fabrie  net,  while  the  net  is  still  stable  with  the 
automatic  position  adjust  of  the  remaining  four  metal  stuffers.  In  Fig.  27-d,  the  final  stability 
stage  is  reached,  and  the  maximum  bending  defleetion  of  the  eomposite  strueture  is  16  mm. 


a)  Stability  stage  A,  all  the  fiber  threads  are  intaet;  the  maximum 
deflection  is  2.5  mm. 
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Broken  thread 


master  thread 


b)  Stability  stage  B,  the  first  master  thread  is  broken,  the  first 
neighboring  thread  beeomes  the  master  thread,  the  maximum 
defleetion  is  4.3  mm. 
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Broken  threads 

e)  Stability  stage  C,  the  first  and  seeond  master  threads  are  broken,  the 
seeond  neighboring  thread  beeomes  the  master  thread,  the  maximum 
defleetion  is  8.3  mm. 


d)  Stability  stage  D,  all  master  threads  are  broken  exeept  the  third 
neighboring  thread  now  beeomes  the  master  thread;  the  maximum 
bending  defleetion  is  reaehed  as  16.0  mm. 

Figure  27.  The  multi-stage  stability  of  the  composite  structure  during  impact. 

The  reaetion  foree  on  the  impaet  objeet  is  shown  in  Fig.  28.  In  the  four  stability  stages,  the 
reaetion  foree  in  stage  A  and  stage  B  are  almost  linear  in  terms  of  the  defleetion.  In  the  last  two 
stability  stages,  the  eomposite  strueture  ean  still  provide  suffieient  bending  stiffness.  Figure  28 
evidenees  the  existenee  of  multi-stage  stability  and  the  effeetiveness  of  the  fiber  ropes  in  the 
eomposite  strueture. 
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Bending  Deflection  (mm) 


Figure  28.  Reaction  force  on  the  impact  object  versus  the  bending  deflection  of  the 

composite  structure 


4.2.5  Future  development 

Figure  29  illustrates  a  complete  fiber-reinforced  composite  model  to  be  developed  in  next  step, 
which  will  be  used  to  predict  material  behaviors  at  the  meso-level.  A  general  homogenization 
method  will  be  then  developed  for  predict  effective  material  properties  of  composite  for  structural 
analysis  and  design  at  the  macro-level.  The  developed  new  capabilities  will  be  finally 
implemented  into  the  FOMD  system  for  optimal  design  of  composite  structures  for  a  wide  range 
of  different  applications. 
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Figure  29.  A  completed  design  model  for  fiber-reinforced  composites 
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4. 3  Potential  Applications 

Figure  30  shows  a  snapshot  of  the  damage  and  deformation  simulation  for  a  HMMWV  under  a 
mine  attaek.  The  landmine  is  detonated  under  the  driver  seat.  As  ean  be  seen  in  Fig.  30,  the  floor 
under  the  seat  has  been  signifieantly  damaged.  The  resultant  aeeeleration  on  the  dummy  ean  also 
been  predieted.  As  in  this  ease,  the  dummy  will  be  seriously  injured  if  no  additional  proteetion  is 
provided. 


Figure  30.  Example  result  of  a  landmine-HMMWV-dummy  blast  simulation. 


Blast  simulations  have  also  been  earried  out  for  investigating  blast  effeets  on  speeial 
substruetures  in  the  vehiele  system.  Figure  31 -a  demonstrates  that  speeial  blast  effeet  -  the 
window  glass  blowout  under  the  blast  wave  ean  be  simulated  using  the  simulation  model.  It  is 
seen  that  onee  the  blast  wave  arrives,  the  glass  is  broken  to  pieees  and  blown  out  from  the 
window,  whieh  matehes  the  aetual  test  results.  Figure  31-b  further  illustrates  that  the  eapability 
eould  be  used  to  design  a  proteetive  armor  between  the  glass  and  driver  for  an  existent  vehiele 
system  against  the  glass  fragments  in  the  blast  event. 


a)  Glass  blowout  under  the  blast  wave 


b)  Glass  fragments  penetrating  an  armor  plate 


Figure  31.  Window  glass  blowout  under  a  blast. 
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For  mine  blast  protection,  structural/material  design  can  be  roughly  divided  into  two 
categories:  1)  design  of  the  protection  panel  (steel/composite)  for  blast  wave  deflection,  and  2) 
material  design  for  blast  energy  absorption. 

In  order  to  protect  the  crew  members  and  the  integrity  of  the  main  structure,  protection 
panels  are  inserted  in  military  vehicles  to  deflect  the  blast  wave.  Figure  32  illustrates  an  example 
design  of  the  centerline  “V-deflector”  assembly,  including  “wing”  deflectors  that  are  added  to 
ensure  adequate  deflection  of  the  blast- wave  away  from  the  cab.  Note  that  a  Crew  Protection  Kit 
(CPK)  prototype  has  been  developed  under  DARPA  support  [21]. 


Crew  Cab 


Figure  32.  Example  design  of  wing/wedge  deflector  geometry  [22]. 

The  successfully  developed  crew/vehicle  protection  kit  can  increase  crew  survivability  of 
tactical  wheeled  vehicles  subject  to  mine  blast.  However,  these  CPK’s  were  based  on  a 
conventional  steel/aluminum  construction,  and  are  close  to  a  half  ton  in  weight.  In  order  to 
reduce  the  weight,  sandwich  structures  have  been  proposed  and  tested  as  an  alternative  [23]. 
Experiments  have  been  carried  out  for  the  composite  panel  under  blast.  It  is  found  the  panel 
cannot  restrain  the  maximum  deflection  as  much  as  the  original  alloy  panel  [23].  Due  to  the 
deflciency  in  inter-laminar  strength  of  the  sandwich  panel’s  faces  and  the  subsequent 
consequence  from  dynamic  blast  loading,  the  dynamic  deflection  of  the  panel  is  not  satisfactory 
compared  to  the  past  monolithic  panel  test  results.  As  studied  by  Yen  and  Jones  [24],  the  failure 
of  the  composite  panel  has  the  following  modes: 

■  Fiber  failure  in  a  blast  loaded  sandwich  panel. 
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■  Delamination  induced  by  in-plane  shock  wave. 

■  Through-the-thickness  compressive  failure  and  spall  fracture. 

It  is  believed  that  the  deficiency  may  be  corrected  with  the  following  improvements: 

□  Orient  some  volumes  of  fibers  or  other  types  of  solids  in  the  through-the-thickness 
direction; 

□  Optimize  ply  layup,  ply  type,  and  matrix  type  to  lower  transverse  tension. 

Figure  33  illustrates  the  capabilities  that  have  been  developed  by  us  in  a  FOMD-Blast 
system  for  predicting  the  pressure  due  to  the  blast  wave  on  a  composite  plate.  In  addition  to  the 
blast  wave,  using  the  new  (under  developing)  mine-soil-composite  interaction  model,  the 
inhomogeneous  property  of  soil,  and  the  interactions  between  the  blast  wave  and  soil  debris  can 
be  accounted  for  in  the  blast  load  prediction  so  that  the  composite  structure  can  be  optimized  with 
more  realistic  loading  conditions.  A  statistical  algorithm  will  be  developed  for  various  soil 
behaviors  under  the  land  explosive,  including  soil  texture,  structure,  failure  properties,  wet/dry 
condition,  etc.  The  effects  of  burial  depth  and  debris-blast  wave  interaction  will  also  be  studied 
to  calibrate  the  dependence  in  blast  energy  capture  of  the  structure. 


Figure  33.  Blast  wave  pressure  modeling  in  the  FOMD-Blast. 
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Figure  34  illustrates  the  FOMD-Blast  main  menu  for  defining  eomposite  layout,  whieh  is 
used  to  seleet  raw  materials  in  the  BTR-Bl  eomposite  and  lay  out  the  eomposition  eonfiguration. 
As  shown  in  Figure  34,  the  BTR-Bl  eomposite  ean  be  eonfigured  with  a  wide  ehoiee  of  different 
raw  materials  for  the  skin,  eore  plate,  stuffer,  fiber,  matrix,  ete.,  and  the  designer  ean  have  great 
freedom  to  ereate  an  innovative  eonfiguration  of  the  BTR-Bl  eomposite  with  needed  effieieney 
and  eonvenienee. 
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Figure  34.  FOMD-Blast  main  menu  for  composite  layout:  selecting  raw  materials  and 
configuration. 


The  present  armor  eoneept  to  defeat  light  to  medium  threats  typieally  eonsists  of  a  hard- 
strong  frontal  surfaee  and  a  softer-tough  and  strong  baeking.  In  addition  to  ballistie  and  blast 
proteetions,  other  major  requirements  for  advaneed  vehiele  armor  are  lightweight,  flexibility, 
maintainability  and  redueed  life-eyele  eost.  Lightweight  is  erueial  to  maintaining  exeellent  road 
and  eross-eountry  mobility,  whieh  is  direetly  related  to  military  deployability  and  survivability. 
Lightweight  is  also  erueial  to  the  transportability  and  sustainability  as  well  as  to  struetural 
integrity  and  durability.  Flexibility  means  the  armor  strueture  ean  be  shaped  or  formed  to  fit 
various  vehiele  eontours.  Maintainability  implies  two  things:  1)  the  integrated  armor  system  ean 
be  easily  installed  and  removed  from  the  vehiele  with  minimal  time  and  manpower,  and  2)  the 
armor  ean  be  easily  repaired  during  war  time  without  replaeing  the  whole  armor  strueture.  Life- 
eyele  eost  is  direetly  related  to  affordability  by  the  US  military  and  the  wide  applieation  of  the 
armor  system.  The  new  integrated  eomposite  design  eapabilities  ean  be  used  to  design  truly 
optimum  and  innovative  eomposite  struetures  for  suiting  the  above  needs. 
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5  CONCLUSIONS 

Meso-structure  of  textile  eomposites  eontributes  to  the  outstanding  properties  of  the  eomposites, 
and  at  the  same  time  beeomes  great  ehallenge  for  simulation  and  design  with  needed  aeeuraey 
and  fidelity.  We  proposed  to  developed  unique  eapabilities  for  modeling  and  simulation  of  textile 
fabries  and  textile  eomposites  based  on  a  truly  nonlinear  model  developed  by  us.  A  damage 
model  for  textile  eomposites  will  be  then  developed  based  on  the  meso-struetural  behaviors  of  the 
eomposite  and  the  use  of  an  enhaneed  homogenization  method.  The  proposed  researeh  will 
eontribute  signifieantly  to  the  related  Army  researeh  areas  for  aehieving  the  Army  next  ten  years 
goals. 
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